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We discuss some aspects of cold atomic gases in the unitarity limit that are of interest 
in connection with the physics of dense hadronic matter. We consider, in particular, the 
equation of state at zero temperature, the magnitude of the pairing gap, and the phase 
diagram at non-zero polarization. 

1. Introduction 

The QCD phase diagram contains several strongly interacting quantum gases or 
liquids. Experiments at the Relativistic Heavy Ion Collider (RHIC) indicate that 
at temperatures close to the critical temperature T c the quark gluon plasma is 
strongly interacting. The strongly interacting plasma (sQGP) is characterized by a 
very small viscosity to entropy ratio and by its large opacity for high energy jets 
1 . In the opposite regime of large baryon density and small temperature the phase 
diagram features several strongly interacting quantum liquids. Electrically neutral 
matter at very low density contains mostly neutrons. Dilute neutron matter has 
positive pressure and is stable even at very low density. However, the neutron- 
neutron scattering length is very large and as a consequence neutron matter is 
strongly correlated even if the density is low 2 . At densities on the order of several 
times nuclear matter density hadronic matter undergoes a phase transition to color 
superconducting quark matter. At extremely large density this phase is weakly 
coupled, but at densities near the transition to nuclear matter the gap and the 
critical temperature are probably large 3 . 

There are many questions about strongly interacting phases of QCD that have 
attracted a lot of interest. Some of these questions are: What are the relevant 
degrees of freedom? Are quasi-particle pictures appropriate? What are the transport 
properties? Is it true that there is a universal bound on the viscosity, and is this 
bound saturated in any of the strongly interacting phases of QCD? 

In this contribution we shall study a simpler system in which many of the 
same questions can be addressed. We consider a cold, dilute gas of fermionic atoms 
in which the scattering length a of the atoms can be controlled experimentally. 
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Fig. 1. The left panel shows the scattering length of 40 K Atoms as a function of the magnetic 
field near a Feshbach resonance, from Regal (2005). The right panel shows the nucleon-nucleon 
scattering length in the 1 So channel as a function of the pion mass. The scatter plot indicates the 
uncertainty due to higher order terms in the chiral effective lagrangian. Figure from Savage and 
Beane (2003). 



These systems can be realized in the laboratory using Feshbach resonances, see 4 
for a review. A small negative scattering length corresponds to a weak attractive 
interaction between the atoms. This case is known as the BCS limit. As the strength 
of the interaction increases the scattering length becomes larger. It diverges at the 
point where a bound state is formed. The point a = oo is called the unitarity limit, 
since the scattering cross section saturates the s-wave unitarity bound a = 4ir/k 2 . 
On the other side of the resonance the scattering length is positive. In the BEC 
limit the interaction is strongly attractive and the fermions form deeply bound 
molecules. 

A dilute gas of fermions in the unitarity limit is a strongly coupled quantum 
liquid that exhibits many interesting properties. On a qualitative level these prop- 
erties are of interest in connection with other strongly interacting field theories. 
The unitarity limit is also quantitatively relevant to the physics of dilute neutron 
matter. The neutron- neutron scattering length is a nn — — 18 fm and the effective 
range is r nn — 2.8 fm. This means that there is a range of densities for which the 
interparticle spacing is large compared to the effective range but small compared 
to the scattering length. It is interesting to note that the neutron scattering length 
depends on the quark masses in a way that is very similar to the dependence of 
atomic scattering lengths on the magnetic field near a Feshbach resonance 5 , see 
Fig. 1. 

2. Universal equation of State 
2.1. Cold Atomic Gases 

We first consider the equation of state of cold fermions in the unitarity limit. We 
are interested in the limit (kpa) — > oo and (fci?r) — > 0, where kp is the Fermi 
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Fig. 2. Lattice results for the energy per particle of a dilute Fermi gas from Lee & Schafer (2005). 
We show the energy per particle in units of 3Ep/5 as a function of temperature in units of Tp. 



momentum, a is the scattering length and r is the effective range. From dimensional 
analysis it is clear that the energy per particle has to be proportional to energy per 
particle of a free Fermi gas at the same density 




where kp is the Fermi momentum. The calculation of the dimcnsionlcss quantity 
£ is a non-perturbative problem. We shall tackle this problem using a combination 
of effective field theory and lattice field theory methods. We first observe that in 
the low density limit the details of the interaction are not important. The physics 
of the unitarity limit is captured by an effective lagrangian of point-like fermions 
interacting via a short-range interaction. The lagrangian is 

£ = ^(^ + £)^T^ 2 ' (2) 

where m is the mass of the fermion and Co is the strength of the four-fcrmion 
interaction. In the weak coupling limit Co is related to the scattering length by 
Co = 47ra/m. 

If the scattering length is small then the energy of the many body system can 
be computed as a perturbative expansion in (ijra). If (fc^a) is large we have to rely 
on numerical simulations. The partition function can be written as 6 



Z = DsDcDc* exp [~S] , 



(3) 



4 Schdfer 



oo 





Fig. 3. Particle-particle ladder diagrams. 

where S is a discrctizcd cuclidean action 

S = ^ [e-^ a *c*(n)c^(fi + 6) - e V-Coa tS (ri)+Sp (1 _ 6/l ) c *(^) c ^)" 

-ft 2 [c*(n)c:(n + f s ) + C *(n)cUn-/ s )]+^ S 2 (n). (4) 

n,l s ,i n 

Here, s is a Hubbard-Stratonovich field, c is a Grassmann field, at is the ratio of 
the temporal and spatial lattice spacings and h = a t /(2m). The sums are over 
spin labels i, lattice sites n, and spatial unit vectors I. is a temporal unit vector. 
Note that for Co < the action is real and standard Monte Carlo simulations are 
possible. Results in the unitarity limit are shown in Fig. 2. From these simulations 
we concluded that £ = (0.09 — 0.42). Lee performed canonical simulations at T = 
and obtained 7 £ = 0.25 . Green Function Monte Carlo calculations give 8 £ = 0.44, 
and finite temperature lattice simulations have been extrapolated to T = to yield 
similar results 9 ' 10 . 

It is also interesting to find analytical approaches to the equation of state in the 
unitarity limit. Since the two-body interaction is large it is natural to begin with 
the sum of all two-body diagrams, see Fig. 3. This sum gives 11 

E = _4 f3 2(fc F a)/(37r) \ 

N 2M \5 + 1-3^(11 -2 Iog(2))(fc F o)J" [> 

from which we deduce £ ~ 0.32. This is reasonably close to the numerical results, 
but since the system is strongly correlated there is no obvious reason to restrict 
ourselves to two-body ladders. We have recently studied the possibility that equ. (5) 
can be justified as the leading term in an expansion in 1/d, where d is the number of 
spatial dimensions 12 > n . This approach appears promising, but 1/d corrections have 
not been worked out yet. In order to obtain a smooth d — > oo limit the coupling 
constant has to be scaled in a specific way. Nussinov & Nussinov observed that if 
the limit a — > oo is taken first then the point d = 4 is special 13 . In d = 4 the 
two-body wave function at a = oo has a 1/r 2 behavior and the normalization is 
concentrated near the origin. This implies that the many body system is equivalent 
to a gas of non-interacting bosons and £(d = 4) = 0. Nishida and Son computed £ 
in d = 3 using an expansion around this limit. At next-to-leading order they find 
14 £ = 0.475. 
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2.2. Strongly Coupled Gauge Theory 

In connection with the RHIC program we are interested in understanding the quark 
gluon plasma in the vicinity of the critical temperature. From lattice simulations 
we know that the energy density quickly reaches about 80% of the ideal gas value, 
and that this ratio is only very weakly temperature dependent on temperature 
for T > 2T C . This behavior can be understood in hard thermal loop rcsummed 
perturbation theory 15 . In this framework the degrees of freedom are dressed quasi- 
quarks and quasi-gluons, and these quasi-particles are weakly interacting. 

Transport properties of the plasma indicate that this may not be the end of the 
story. If quasi-particles are indeed weakly interacting then it is hard to see how the 
viscosity can be anomalously low or the opacity be very large. A complementary, 
strong coupling, calculation was performed in the strong coupling and large N c 
limit of Af = 4 SUSY Yang Mills theory. The calculation is based on the duality 
between the strongly coupled gauge theory and weakly coupled string theory on 
AdSs x £5 discovered by Maldacena 16 . The correspondence can be extended to 
finite temperature. In this case the relevant configurations is an AdS$ black hole. 
The temperature of the gauge theory is given by the Hawking temperature of the 
black hole, and the entropy is given by the Hawking-Beckenstein formula. The result 
is that the entropy density of the strongly coupled field theory is equal to 3/4 of the 
free field theory value 17 . Clearly, the number 3/4 can be viewed as a gauge theory 
analog of the parameter £ studied in the previous section. The remarkable result 
is that the gauge theory value is so close to one, so that based on the equation of 
state alone one cannot easily distinguish a strongly and a weakly coupled system. 

3. Pairing 

3.1. Cold Atomic Gases 

At low temperature the atomic gas becomes superfluid. If the coupling is weak then 
the gap and the critical temperature can be calculated using BCS theory. The result 
is 



where the factor (4e) 1//3 is the screening correction first computed by Gorkov et 
al. 18 . Higher order corrections are suppressed by powers of (fcffl). In BCS theory 
the critical temperature is given by T c — e^A/ir. Clearly, the critical temperature 
grows with the scattering length. Naively extrapolating equ. (6) to the unitarity 
limit gives T c ~ 0.28-Ep- In the BEC limit the critical temperature is given by the 
Einstein result T c = 3.31/?^ /tom, where pu — p/2 and tom = 2m are the density 
and mass of the molecules. This implies that T c ~ 0.21Ep where we have defined the 
Fermi energy as Ep = k F /(2m) = (3tt 2 p) 2 ^ 3 / (2m). Interactions between the bosons 
increase the critical temperature 19 , T c = T° + 0(aBP 1 ^ 3 )- Near the unitarity limit 
the boson scattering length as is proportional to the scattering length between the 




(6) 
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fermions 20 , cibb — 0.6a. These results suggest that the unitarity limit corresponds 
to the maximum possible value of T C /E F . 

The value of T c has been determined in a number of lattice calculations. A careful 
scaling analysis by Burovski et al. yields 10 T c = 0.152(7)Ep. Earlier determinations 
of T c can be found in 21,9 . A direct calculation of the pairing gap at zero temperature 
using Green Function Monte Carlo methods gives 22 A = Q.bAEp. 

3.2. Color Superconductivity 

Pairing also takes place in the high density, low temperature phase of QCD. At 
asymptotically large density the attraction is due to one-gluon exchange between 
fermions with opposite momenta and anti-symmetric spin, color and flavor wave 
functions. In this limit the pairing gap is given by 23 ' 24 

A = 2A BCS exp(-^±i)e X p(-^). (7) 

where g is the running coupling constant evaluated at the scale /x and Abcs — 
256tt 4 (2/N /) 5 / 2 g~ 5 u. Here, u is the baryon chemical potential and Np is the num- 
ber of flavors. This result exhibits a non-BCS like dependence on the coupling 
constant which is related to the presence of unscreened magnetic gluon exchanges. 
The critical temperature is nevertheless given by the BCS result T c = e^A/ir. 

If the weak coupling limit the gap and the critical temperature are exponentially 
small. The ratio T c /Ep increases with g and reaches a maximum of T c = 0.025Ep 
at g — 4.2. The maximum occurs at strong coupling and the result is not reliable. 
Using phcnomcnological interactions, or extrapolating the QCD Dyson-Schwinger 
equations into the strong coupling domain 25 , one finds critical temperatures as large 
as T c — 0.15Ep. In connection with phcnomcnological applications to neutron stars 
and the physics of heavy ion collisions in the regime of the highest baryon densities 
it is very important to reduce the uncertainty in these estimates. One possibility 
is to use lattice studies of QCD-like theories in which the euclidean action remains 
positive at finite chemical potential. Examples arc QCD with two colors and QCD 
at finite isospin density. 

4. Stressed Pairing 

4.1. Polarized Cold Atomic Gases 

The superfluid state discussed in Sect. 3.1 involves equal numbers of spin up and 
spin down fermions. One aspect of the paired state that has attracted a lot of 
interest is the response of this system to a non-zero chemical potential coupled 
to the third component of spin, 5u = u^ — u^. A conjectured (and, most likely, 
oversimplified) phase diagram for a polarized gas is shown in Fig. 4. In the BEC 
limit the gas consists of tightly bound spin singlet molecules. Adding an extra up 
or down spin requires energy A. For \8jj\ > A the system is a homogeneous mixture 
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Fig. 4. Conjectured phase diagram for a polarized cold atomic Fermi gas as a function of the 
scattering length a and the difference in the chemical potentials 5/J, = /i* — fj,\ , from Son & 
Stephanov (2005). 

of a Bose condensate and a fully polarized Fermi gas. One can show that in the 
dilute limit this mixture is stable with regard to phase separation 2 . 

We can also analyze the system in the BCS limit. This analysis goes back to the 
work of Larkin, Ovchninikov, Fulde and Ferell (LOFF) 27 > 28 J see the review 29 . We 
first consider homogeneous solutions to the BCS gap equation for <5/j. 7^ 0. In the 
regime Sfi < Ao where Ao = A(5/i = 0) the gap equation has a solution with gap 
parameter A = A . This solution is stable if Sfj. < A a /^/2 but only meta-stable in 
the regime Aq/^/2 < Sfi < Ao. The BCS solution has vanishing polarization. The 
transition to a polarized normal phase is first order, and systems at intermediate 
polarization correspond to mixed phases. 

LOFF studied whether it is possible to find a stable solution in which the gap 
has a spatially varying phase 



This solution exists in the LOFF window Sfii < <5/i < <5/i2 with S/ii = Ao/y/2 ~ 
0.7lAo and <5/i2 — 0.754Ao. The LOFF momentum q depends on 5fi. Near 6/12 we 
have qvF — 1.25/i, where vf is the Fermi velocity. The gap A goes to zero near 8^2 
and reaches A ~ 0.25Ao at S(i±. 

These results suggest that at some point on the phase diagram between the BEC 
and BCS limits the homogeneous superfluid becomes unstable with respect to the 
formation of a non-zero supercurrent Vip, where <p is the phase of the condensate. 
We can study this question using the effective lagrangian 



A(x) = Ae 2l ^ s . 



(8) 



C = (id, - e(-W) - i0<p) - A) ^ + h. - L 0^ 



(9) 



8 Schdfer 



Here, tjj describes a fermion with dispersion law e{p) and <p is the superfluid Gold- 
stone mode. The low energy parameters f t and / are related to the density and the 
velocity of sound. The p-wave coupling of the fermions to the Goldstone boson is 
governed by the U(l) symmetry of the theory 

Setting up a current v s — dtp/m requires energy f 2 m 2 v 2 /2. The contribution 
from fermions can be computed using the fermion dispersion law in the presence of 
a non-zero current 

tv{p) = e(p) + v s -p- Su . (10) 

The total free energy is 

F(v s ) = \nmvl + j ^ e v (p)Q (-e v (p)) , (11) 

where n is the density and we have used / = n/m. Son and Stephanov noticed 
that the stability of the homogeneous phase depends crucially on the nature of 
the dispersion law e(p). For small momenta we can write e(p) ~ e + a P 2 + ftp 4 - 
In the BEC limit a > and the minimum of the dispersion curve is at p = 
while in the BCS limit a < and the minimum is at p ^ 0. In the latter case 
the density of states on the Fermi surface is finite and the system is unstable with 
respect to the formation of a non-zero current. On the other hand, if the minimum 
of the dispersion curve is at zero, then the density of states vanishes and there is 
no instability. As a consequence there is a critical point along the BEC-BCS line 
at which the instability will set in. 

Clearly, the LOFF solution is of the same type as the supercurrent state. The 
difference is that in the supercurrent state the current is much smaller than the 
gap, vp(\7(p) <c A, while in the LOFF phase vf (V<^) > A (and vp(\7(p) ^> A near 
S^2 )• In the supercurrent state the Fermi surface is mostly gapped but a small shell 
near one of the pole caps is ungapped. In the weakly coupled LOFF state there are 
gapless excitations over most of the Fermi surface but pairing takes place near two 
rings on the northern and southern hemisphere. 

4.2. CFL Phase at Non-zero Strange Quark Mass 

The exact nature of the color superconducting phase in QCD depends on the baryon 
chemical potential, the number of quark flavors and on their masses. If the baryon 
chemical is much larger than the quark masses then the ground state of QCD with 
three flavors is the color-flavor-locked (CFL) phase. The CFL phase is characterized 
by the pair condensate 30 

WC"y&$> = - (12) 

This condensate leads to a gap in the excitation spectrum of all fermions and 
completely screens the gluonic interaction. Both the chiral SU(3)l x SU(3)r and 
color SU (3) symmetry are broken, but a vector-like SU(3) flavor symmetry remains 
unbroken. 
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In the real world the strange quark mass is not equal to the masses of the up 
and down quark and flavor symmetry is broken. At high baryon density the effect 
of the quark masses is governed by the shift /j, q = m^/(2ji) of the Fermi energy 
due to the quark mass. This implies that flavor symmetry breaking becomes more 
important as the density is lowered. Clearly, the most important effect is due to the 
strange quark mass. There are two scales that are important in relation to fi s . The 
first is the mass of the lightest strange Goldstone boson, the kaon, and the second 
is the gap A for fermionic excitations. When fi s becomes equal to tuk the CFL 
phase undergoes a transition to a phase with kaon condensation 31 . In the following 
we shall study the phase structure near fi s ~ A. 

Our starting point is the effective theory of the CFL phase derived in 32 > 33 . The 
effective lagrangian contains Goldstone boson fields £ and baryon fields N. The 
meson fields arise from chiral symmetry breaking in the CFL phase. The leading 
terms in the effective theory are 

C = ^Tr (VoSVoSt - t^VEVEt) , (13) 

where is the pion decay constant. The chiral field transforms as S — > LY*R) 
under chiral transformations (L,R) e SU(3)l x SU(3)r. Baryon fields originate 
from quark- hadron complementarity 34 . The effective lagrangian is 

C = Tr (NUv^D^N) - DTt (NV75 {A^ N}) - FTr (NV75 [A^,N]) 

+ j { (Tr (N L N L ) - [Tr (A7J] 2 ) - (Tr (N R N R ) - [Tr (A^)f) + h.c} , (14) 

where Nl^r are left and right handed baryon fields in the adjoint representation 
of flavor SU(3), = (1, v) is the Fermi velocity, and A is the superfluid gap. We 
can think of N as being composed of a quark and a diquark field, Nl ~ 9l(<?l9l)- 
The interaction of the baryon field with the Goldstone bosons is dictated by chiral 
symmetry. The covariant derivative is given by D^N = d^N + «[V M , N}. The vector 
and axial-vector currents are 

V M = ~{^ t + ^}, A, = - l -Z(d^)H, (15) 

where £ is defined by £ 2 = X. The low energy constants f^,v v ,D,F can be cal- 
culated in perturbative QCD. Symmetry arguments can be used to determine the 
leading mass terms in the effective lagrangian. Bedaque and Schafer observed that 
Xl = MM^ /{2pp) and X R = M^M/(2pp) act as effective chemical potentials and 
enter the theory like the temporal components of left and right handed flavor gauge 
fields 31 . We can make the effective lagrangian invariant under this symmetry by 
introducing the covariant derivatives 

D N = d Q N + i[T , N], T = - % - {£ (do + iX R ) g + ? (d + iX L ) £} , (16) 
V £ = d £ + iX L £-z£A^. (17) 
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Fig. 5. Left panel: Energy density as a function of the current jk for several different values 
of fi s = m1/(2pp) close to the phase transition. Right panel: Ground state energy density as 
a function of fi s . We show the CFL phase, the kaon condensed CFL (KCFL) phase, and the 
supcrcurrent state (cur KCFL). 



Using equ. (14-16) we can calculate the dependence of the gap in the fermion 
spectrum on the strange quark mass. For m s — there are 8 quasi-particles with 
gap A and one quasi-particle with gap 2A. As m s increases some of the gaps 
decrease. In the K° condensed phase the gap of the lowest mode is approximately 
given by A = Ao — 3/i s /4 where jj, a — m 2 s l{2pp) and Ao is the gap in the chiral 
limit. 

For fi s ~ 4Ao/3 the system contains low energy fermions derivatively coupled 
to Goldstone bosons. The situation is essentially equivalent to the cold atomic 
system studied in the previous section. The main difference is that in the CFL 
phase there are many more currents that can appear. In the K° condensed phase 
the natural ansatz is a hypercharge current carried by the kaon field. We take 
£(x) — U{x)^ K oW{x) where ^ K o is the K° background, U(x) is a hypercharge 
transformation and jk = iU^WU is the kaon current. The dispersion relation of the 
lowest mode is 35,36 

I 2 3 1 

^ = Ao+ 2Ao" _ 4 Ms ~4 iT ' Jk ' (18) 

where I is the momentum relative to the Fermi surface. The energy relative to the 
CFL phase is the kinetic energy of the current plus the energy of occupied gapless 
modes 

£ = \vlfJl + '^JdiJ§ unO(-un). (19) 

The energy functional can develop a minimum at non-zero jk because the current 
lowers the energy of the fermions near one of the pole caps on the Fermi surface. 
Introducing the dimensionless variables x — jk/ (o,A) and h — (3/i s — 4A) / (aA) we 
can write 

£ = Cf h (x), f h (x) =x 2 -- \(h + xf' 2 Q{h + x) -(h- xf/ 2 Q(h - x)] , (20) 

X L J 



From Trapped Atoms to Liberated Quarks 11 



where C and a are numerical constants. The functional given in equ. (20) was 
analyzed in 37 ' 35 > 36 j se e Fig. 5. There is a critical chemical potential fi s = (4/3 + 
ah cr it/3)A above which the groundstate contains a non-zero supercurrent jk- This 
current is canceled by a backflow of gapless fermions. The Goldstone current phase 
is analogous to the supercurrent state in cold atomic systems, and to p-wave pion 
condensates in nuclear matter. 

The Goldstone current phase in pure CFL matter, without an s-wave kaon 
condensate, was analyzed by Gerhold and Schafer 38 . In that case the structure of 
the current is different, but the energy functional and the nature of the instability 
are very similar. Gerhold and Schafer showed, in particular, that the magnetic 
screening masses in the Goldstone boson current phase are real. Near the onset of 
the instability the current is small compared to the gap, but the current grows with 
the mismatch /j, s . Once the current becomes comparable to the gap it may be more 
appropriate to characterize the system as a LOFF state 39 . 

5. Outlook 

There are many questions about cold fermion gases in the unitarity limit, and about 
their relevance to strongly interacting hadronic matter, that remain to be resolved. 
The unitarity limit is most directly connected with the physics of dilute neutron 
matter. In neutron matter the scattering length is not exactly infinite, and the 
effective range is not zero. It is clearly important to understand how corrections 
due to the finite scattering length and effective range affect the equation of state 
and other observables. 

There is a great deal of experimental effort devoted to the study of polarized 
fermionic gases. Mixed phases of superfluids and polarized normal fluids have been 
observed 40 , but so far none of the predicted inhomogeneous or gapless superfluids 
have been experimentally detected. Finally, much work remains to done with regard 
to the transport properties of cold fermionic gases. We would like to understand, 
both theoretically and experimentally, how the shear viscosity depends on the tem- 
perature and the scattering length. Damping of collective oscillations of trapped 
fermions has been observed experimentally 41 , but the damping coefficients have 
not been related to transport properties of the bulk system. 
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